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Transit timing variation in exoplanet WASP-3b 
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ABSTRACT 

Photometric follow-ups of transiting exoplanets may lead to discoveries of additional, 
less massive bodies in extrasolar systems. This is possible by detecting and then 
analysing variations in transit timing of transiting exoplanets. We present photometric 
observations gathered in 2009 and 2010 for exoplanet WASP-3b during the dedicated 
transit-timing-variation campaign. The observed transit timing cannot be explained 
by a constant period but by a periodic variation in the observations minus calcula- 
tions diagram. Simplified models assuming the existence of a perturbing planet in the 
system and reproducing the observed variations of timing residuals were identified by 
three-body simulations. We found that the configuration with the hypothetical second 
planet of the mass of ~15 M^, located close to the outer 2:1 mean motion resonance is 
the most likely scenario reproducing observed transit timing. We emphasize, however, 
that more observations are required to constrain better the parameters of the hypo- 
thetical second planet in WASP-3 system. For final interpretation not only transit 
timing but also photometric observations of the transit of the predicted second planet 
and the high precision radial-velocity data are needed. 
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1 INTRODUCTION 

The analysis of observations minus calculations (O — C) 
diagrams is a great tool commonly used for astrophysi- 
cal studies of eclipsing binaries or pulsating variables. This 
method has already brought discoveries of exoplanets. Vari- 
ations in the O — C diagram of the pulsating subdwarf B 
HS 2201+2 610 were interprete d as a fingerprint of a gi- 
ant planet (|Silvotti et al.l I2007T ) . More recently, two plan- 
ets orbiting an e clips ing system HW V ir were announced by 
iLee et all (j2009h and lQian et"aH (|2010l ) found a giant planet 
around the hibernating cataclysmic binary system QS Vir. 
These discoveries are based on interpretation of the light- 
travel time effect generated by an additional body. 

The timing of known transiting planets is also 
expected to lead to discovering additional planets 
<|Miralda-Escudell2002l : iHolman fc Murravll200 j ; lAgol et all 
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Friedrich Schiller University, 
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120051 : ISteffen et al.|[2007h . In a single-planet extrasolar sys- 
tem a planet orbits its host star on a Keplerian orbit. If 
one assumes that the inclination of its orbit plane is close 
to 90°, transits occur at equal intervals. If there is another 
(not necessarily transiting) planet in the system it inter- 
acts gravitationally with the transiting planet and generates 
deviations from the strictly Keplerian case. These pertur- 
bations result in a quasi-periodic signal in an O — C dia- 
gram of the transiting planet. The transit timing variation 
(TTV) method can be sensitive to small perturbing masses 
in orbits near low -order mean-motion resonances (MMRs, 
ISteffen et alj|2007h . A terrestrial-mass planet perturbing a 
hot- Jupiter gas giant is expected to cause a TTV ampli- 
tude of ~1 minute an d this signal grows sharply as bod- 
ies approach a MMR ( Steffen et al.l 1213071 ). One must note 
that deriving the orbital elements and mass of the per- 
turb er from TTV, however, is a difficult inverse problem 
(e.g. iNesvornv fc Morbidell1 l2008). For a given transiting 
planet the TTV signal depends on 7 unknown parameters 
of the perturbing planet, i.e. its mass, semi-major axis, ec- 
centricity, inclination, nodal and periapse longitudes and dif- 
ference of mean orbital phases of both planets for a given 
epoch. Exploring su ch a 7-dimensional space o f parameters 
is not a trivial task (Ne svornv fc B eauge 20ld). The differ- 
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Table 1. The summary of observing runs: X - the airmass range, iV eX p — the number of 
useful exposures, T cxp — used exposure times, FWHM - the averaged full width at half 
maximum of a stellar profile. Dates are given at the beginning of nights. 



Run 


Date 


Observatory 


X 




N exp 




FWHM (") 


1 


2009 July 28 


Rozhen 


1.01 - 


1.29 


283 


30 


1.5 


2 


2009 Aug. 21 


Rozhcn 


1.04 - 


1.92 


200 


60 


1.5 


3 


2009 Sep. 03 


Rozhen 


1.01 - 


1.38 


298 


30 


1.3 


4 


2009 Sep. 27 


Jena 


1.05 - 


1.94 


797 


8 


3.2 


5 


2009 Nov. 03 


Jena 


1.13 - 


2.32 


531 


10, 15 


3.9 


6 


2010 Apr. 18 


Jena 


1.08 - 


2.19 


412 


25, 30 


4.8 



ent configurations may generate si milar TTV signals wi th 
the identical dominant periodicity (|Ford fc Holman 2007). 

The timing effects have been alr eady studied in sev- 
eral transiting pla net system: Tr ES-1 JSteffen fc Ago]||2005| : 
Rabus et all 120091), H P 189733 (|Miller-Ricci et all I2008b1 
Hrudkova et~ai1 l2010h. H P 209458 (lAgol fc Steffenl 120071 



rirudKova et al.l izuiui). 1 11J i!uy40o_{AKOl &c stcttcn zvvt; 
Miller-Ricci et al.ll2008al ), GJ 4 36 dBean fc Seifahrdl200^T . 



C0R0T-I (lBeanll2009h . TrES-2 dRabus et al.ll2009l), TrE S-3 
l|Gibson et al.ll2009h . and HAT-P-3 l|Gibson et al.ll2010l ). In 

all these cases the authors could only put constraints on the 
parameters of a hypothetical second planet because no sig- 
nificant signals in O — C diagrams have been found so far. 



iPiaz et alT (|2008l ) announced detection of period variation 
in OGLE-TR-lllb and suggested that it may be caused by 
a perturbing Earth- mass planet in an outer orbit. However, 
lAdams et all l|2010t ) observed 6 new transits of OGLE-TR- 
lllb and concluded that there is no compelling evidence for 
the timing variation in the system and placed an upper limit 
of about 1 M0 on the mass of a potential second planet in 
a 2:1 orbital resonance. 

Our first efforts to detect TTV signal for the transiting 
planets XO-lb, TrES-lb and TrE S-2b resulted in re deter- 
mining their transit ephemerides (|Raetz et al.ll2009al fbTl. In 
2009 we launched an international observing campaign to 
detect and characterise a TTV signal in selected transiting 
exoplanets. The programme is realised by collecting data 
from 0.6-2.2-m telescopes spread worldwide at different lon- 
gitudes. In this paper we present photometric observations 
obtained for transiting exoplanet WASP-3b. 

The existence of a planet around the unevolved main 
sequence star WASP-3 (V = 10.485 mag) was discovered 
by det ecting transits follow ed by radial velocity measure- 
ments |Pollacco et al.l [20081 ) . The exoplanet turned out to 
be a strongly irradiated gas giant. Its mass and radius were 
found to be 1.76±£?f Mj and 1.3lt^;°4 Rj, respectively. It 
orbits its host star in 



^44 h in a circular orbit with the 



semi-major axis of 0.0317 



+0.0005 



au. The orbital inclination 



was found to be 84.4^' g deg. The host star has a spectral 
ty pe of F7-8V and pho tospheric temperature of 6400 ± 100 
K. iGibson et al.l ()2008l ) observed two transits with a photo- 
metric precision of ~4 millimag (mmag) and redetermined 
the planetary radius 1.291q'i2 an d the orbital inclina- 
tion 85.06t ;is deg. The mean d ensity of the planet was 
found to be 0.821°;^ pj. Moreover. IGibson et al] (120081 ) cal- 
culated new ephemeris which resulted in the orbital period 
of 1 .846835 ± 0.000002 d 

ISimpson et al.l l|2010T ) presented a spectroscopic obser- 
vation of the Rossiter-McLaughlin effect for WASP-3 sys- 
tem. The sky-projected angle between the stellar rotation 



axis and planetary orbital axis A was found to be 15lg° 
deg, i.e. consi s tent with zero within 2a. More recently, 
ITripathi et all (|2010h obtained spectra during two sepa- 
rate transits and got a more precise value A = 3.3I44 
deg. These findings indicate that WASP-3b underwent rela- 
tively non-violent migration process which did not perturb 
it from the primordia l alignment of th e proto-planeta r y disk 
l|Simpson et al.ll2010l ). Furthermore, ITripathi et"ai] feoich 
presented photometric observations of 6 transits, including 
one incomplete. Some deviations from the linear ephemeris 
were found and explained as a results of either genuine pe- 
riod vari ations or unde r estim ating the transit time uncer- 
tainties. ITripathi et al.l (|2010l ) also redetermined planets's 
mass which was found to be 2.04 ± 0.07 Mj. 

The order of this paper is as follows. In Section 2 ob- 
servations and data reduction are presented. The finding 
that the observed transit timing cannot be explained by a 
constant period is shown in Section 3. Preliminary configu- 
rations including a third body in the system (a perturbing 
planet), which may reproduce the observed variations of tim- 
ing residuals, are identified and discussed in Section 4. Final 
conclusions are collected in Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed 6 complete transits of WASP-3b with two tele- 
scopes both with the effective mirror diameter being 60 
cm. While preparing observing runs only predicted com- 
plete transits during darkness were considered. Observations 
started ~1 hour before the expected beginning of a tran- 
sit and ended ~1 hour after the event. Collecting long se- 
ries of out-of-transit data is crucial in r emoving systematic 
trends and normalizi ng a light curve (e.g. IGibson et ai1l2009l ; 
iHrudkova et aL 201C ). All observations were acquired in the 
.fi-band filter (|BessellTl99Ct ). 



2.1 Rozhen 

The CCP photometry of 3 transits on 2009 July 28, August 
21 and September 3 was gathered with the 60-cm Cassegrain 
telescope at the Rozhen National Astronomical Observatory 
(NAO, Bulgaria), operated by the Institute of Astronomy, 
Bulgarian Academy of Sciences. The telescope was equipped 
with the CCP camera FLI PL09000 (3056 x 3056, 12pm 
pixel). The field of view was 17!3 x 17^3 with the scale of 
0.334 arcsec per pixel. As a result of 3 x 3 binning the final 
scale was 1.017 arcsec per pixel. Exposure times were se- 
lected on nightly basis to maximize the signal-to-noise (S /N) 
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Figure 1. Light curves of WASP-3b transits in individual runs. The best-fitting models are shown as continuous lines. The residuals are 
presented in panels below each light curve. 



ratio, keeping pixels unsaturated and intensity in the linear 
region of the CCD (<50000 counts) for the object and prob- 
able standard stars. 

The standard IDL procedures (adapted from 
DAOPHOT) were used for the reduction of the pho- 
tometric data (dark frame subtraction and flat- fielding) , 
and computing the differential aperture photometry. To 
measure instrumental magnitudes, we used apertures 
with radii of 3.2, 4, 5, 6, 7, 8 and 9 binned pixels and a 
background ring located between 12 and 16 pixels. The best 
precision was reache d using an apertur e of 4 binned pixels. 
Using the method of lEverettT fc Howell (200l|) several stars 
(4 to 6) with photometric precision better than 5 mmag 
were selected to create an artificial standard star used for 
differential photometry. 



2.2 University Observatory Jena 

Transits on 2009 September 27, November 3 and 2010 April 
18 were monitored with the 90-cm reflector telescope at the 
University Observatory Jena. Observations were taken with 
the new CCD-imager STK (see iMugrauer fe Bertholdl201(il ) 
which was installed in the Schmidt-focus of the telescope 
(D = 60 cm, f/D = 3). The camera exhibits 52'8 x 52'.8 field 
of view with a pixel-scale of 1.546 arcsec per pixel. The expo- 
sure time was set to get the highest S /N ratio for the target 
and was adjusted during a run due to changes of the airmass 
and atmospheric transparency. During third run the tele- 
scope was significantly defocused and hence longer exposure 
times were used. Applying this method wa s expected to min- 
imise random and flat-fielding errors (e.g. ISouthworth et all 
2009), as confirmed by our higher-quality light curve (see 
bottom right panel in Fig. [l]). 



CCD frames were processed using a standard proce- 
dure that included subtraction of a dark frame, flat-fielding 
with twilight flats, differential aperture photometry, and 
astrometric calibration with the software pipeline devel- 
oped for the Semi-Automatic Variabilit y Search sky survey 
(Nicdziclski, Macicicwsk i fe Czartl [20031 ) . The instrumental 
coordinates of stars in frames were transformed into equa- 
torial ones making use of positions of stars brighter than 
V — 15.5 mag and extracted from the Guide Star Cata- 
logue. To measure instrumental magnitudes, we used aper- 
tures with radii of 5, 6, 7 and 8 pixels and a background ring 
located between 12 and 20 pixels. The 6 pixel aperture was 
found to produce light curves with the smallest scatter and 
they were used in further analysis. To generate an artificial 
comparison star, 30 per cent of the stars with the lowest 
light-curve scatter were selected iteratively from field stars 
not fainter than 3 magnitudes below the saturation level. 
To achieve the similar probing rate in all runs (~1 point per 
minute), 3 and 2 point binning of a final light curve was ap- 
plied. The detailed summary of observing runs is presented 
in Table Q] 



3 RESULTS 

3.1 Light curve analysis 

A model-fittin g algorithm available via the Ex oplanet Tran- 
sit Database (|Poddanv. Brat fe Peichal |2010T ) was used to 
derive transit parameters: transit duration, depth and mid- 
transit time and their err ors. The proced u re em ploys the 
OCCULTSMALL routine of iMandel fe Agoll l|2002f ) and the 
Levenberg-Marquardt non-linear least squares fitting algo- 
rithm. The latter routine also provides errors for the fitted 
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Table 2. Parameters of transit light-curve modelling. To, Tj, <5, cr, and E denote the mid-transit time, the 
transit time duration, the depth, the averaged standard deviation of the fit, and the epoch, respectively. 
The O — C values are given both in days and in errors of mid-transit times. Note that BJD times are 



based 


on Terrestrial Dynamic 


Time (TT). 












Run 


To 


T d 


5 


a 


E 


O-C 


O-C 




BJD - 2455000 


(min) 


(mmag) 


(mmag) 




(d) 


(To errors) 


1 


41.41271 ± 0.00049 


161.0 ± 1.6 


12.4 ±0.3 


1.6 


236 


-0.00046 


-0.9 


2 


65.41995 ± 0.00059 


158.9 ± 1.9 


12.6 ±0.4 


1.9 


249 


-0.00207 


-3.5 


3 


78.34873 ± 0.00058 


163.7 ± 1.9 


13.9 ±0.5 


2.4 


256 


-0.00114 


-2.0 


4 


102.35933 ± 0.00056 


157.9 ± 1.9 


12.3 ±0.4 


2.1 


269 


±0.00120 


±2.1 


5 


139.29713 ± 0.00049 


162.9 ± 1.6 


13.4 ±0.4 


1.7 


289 


±0.00169 


±3.4 


6 


305.51082 ± 0.00039 


162.6 ± 1.3 


12.1 ±0.6 


1.2 


379 


±0.00018 


±0.5 



parameters (O. Pejcha 2010, private communication). Our 
tests showed that the errors of mid-transit times derived 
from Levenberg-Marquardt method are consistent within 
±30 per cent with these values obtained from a Markov 
Chain Monte Carlo algorithm. An impact parameter b — 
acosi/R* — 0.448 ± 0.014, where a is semi-major axis, i 
is inclination , and 7?* is host-star radius, was taken from 
ICibson et all (2008) and was fixed during the fitting pro- 
cedure. The model-fitting algorithm uses the linear limb- 
darkening law, thus we used a linear limb-darkening law by 
IVan Hammel jl993l ) with the linear limb-darkening coeffi- 
cient linearly interpolated for the host star. To determine 
the zero-point shift of magnitudes and to remove system- 
atic trends which may exist in our data, a first- or second- 
order polynomial was used. The mid-transit times were cor- 
rected from UTC to Terrestrial Dynamic Time (TT) and 
then transformed into BJD. Light curves with best-fitting 
models and residuals are shown in Fig. [T] Derived parame- 
ters are collected in Table [2] 

We achieved an averaged photometric precision between 
1.2 and 2.4 mmag. The mid-transit timing errors are in the 
range of 34-51 s. The mean transit duration was found to be 
161.2±2.3 min - a value similar to 159. Sj^ g min reported by 



iPollacco et al.l (120081') and not i ceably smaller than 165.2+Qg 
min reported by Gibson ct al. (2008) and 168.8±0.7 min an- 
nounced by iTripathi et all |2010l ). The mean transit depth 
was found to be 12.8 ± 0.7 mmag which results in a planet- 
to-star radii ratio p = 0.108 ± 0.003. This value is within 
the range of planet- to-star radii ra ti os de te rmined for in- 
dividu al t ransits by ITripathi et all (|2010l ). IPollacco et al.l 
HH) and lGibson et al.l (|2008l ) obtained smaller values, i.e. 



0.1030 



+0.0010 



and p = 0.1014+2'ooog, respectively. 



3.2 Transit ephemeris 

Having a long time span of observations we determined a 
new ephemeris. As a result of fitting a linear function of 
epoch and period P, we obtained: 

To = 2454605.56000 ± 0.00011 (BJD, based on TT) 

P = 1.8468355 ± 0.0000007 d. 

Individual mid-transit errors were taken as weights. 

The O — C diagram was generated using the new 
ephemeris. The timing residuals are collected in Table[2] The 
O—C diagram is plotted in Fig.[2]where besides transit times 
reported in this paper the literature data are also shown. Fol- 
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Figure 2. The observed minus calculated diagram for WASP- 
3b. The open symb ols denote literature d ata (complete tran - 
sits only) taken from IPollacco et al.l j2008h . iGibson et alj (l2008r i 
and ITripathi et al. while the filled ones denote results re- 

ported in this paper. Timing residuals were calculated using a 
new ephemeris derived in this study. 



lowing lAdams et al] (|2010|) we hav e confirmed that the t imes 
published bylPollacco et~atl (|2008l ), lGibson et all (|2008h and 
ITripathi et al] (|2010h do not account for the UTC-TT cor- 
rection (D. Pollacco 2010, private communication; A. Tri- 
pathi, private communication). Therefore, we added appro- 
priate corrections to published times before analysing the 
O — C diagram. In the case of 3 points (1 from ITripathi et al] 
2010 and 2 reported in this paper) timing residuals have a 
significance greater than 3 a. That indicates that the uncer- 
tainties are underestimated or that the orbital period is not 
constant. 

It is worth mentioning that data points from Rozhen lie 
below or close to zero level in the O — C diagram while those 
from Jena lie above it. This finding could suggest the pres- 
ence of the systematic offset between both telescopes. How- 
ever, our experience shows that this scenario seems to be 
unlikely. Simultaneous observations collected by both tele- 
scopes for another transiting planet gave difference between 
mid-transit times of 10 s only - much less than the error 
bars (Maciejewski et al. 2010, in prep.). One must note that 
error underestimating cannot be completely ruled out in our 
determinations but it would have a negligible influence on 
data point distribution and final conclusions. 
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Figure 3. (a) A periodogram generated for timing residuals plot- 
ted in FigEl showing a peak at /ttv = 0.0145 cycl P -1 . (b) The 
O — C diagram phased with Pttv = 69.0 ± 2.4 P corresponding to 
/ttv- The open and filled symbols denote literature and this pa- 
per (see Table [2]l data, respectively. Note that data points were 
duplicated to cover phases from —0.5 to 1.5. 

3.3 Search for a cyclic TTV 

The timing residuals were searched for cyclic transit time 
variation by means of the analysis of variance method 
(ANOVA. ISchwarzenberg-Czernvl ll996) to identify any sig- 
nificant periodicities. Fig. [3^ presents the generated pe- 
riodogram covering frequencies / in the range of 0.000- 
0.038 cycl P" 1 where the upper limit is the calculated 
Nyquist frequency /n- A significant peak was found around 
/ttv = 0.0145 cycl P _1 which corresponds to a period of 
Pttv = 69.0 ± 2.4 P. The O - C diagram phase-folded at 
Pttv is plotted in Fig. [3p where the best-fitting sinusoid is 
also sketched. The semi-amplitude of variation was found to 
be 0.0014 ± 0.0002 d. 



4 DISCUSSION 
4.1 Exomoon 

A sinusoidal TTV might be caused by an exomoon 
around the transiti n g planet (ISimon. Szatmarv fc Szabdl 
120071 : iKippind [2009h . Kippind (§009) predicts that an ex- 
omoon also produces transit duration variation (TDV) with 
the same period as TTV and shifted by 7t/2 in phase. That 
effect is not visible in TDV of WASP-3b. If the hypothetical 
moon is located at a Hill radius away from the planet, its 
mass is expected to be 6 per cent of the planetary mass to 
produce the observed TTV amplitude. This extreme value is 



a lower limit of the mass because the needed mass increases 
as the moon is closer to the planet. Both tests indicate that 
the exomoon hypothesis is unlikely here. 

4.2 Reanalysis of the radial velocity data 

The possible non-zero eccentricity of the transiting planet 
may be a results of confusion with a two-planet system in 
which both planets orbit its h ost star in circular orbits an d 
stay in an inner 2:1 resonance |Anglada-Escude et al.ll201ul ). 
This degeneracy is a consequence of the Fourier expansion of 
the Kepler equation into powers of eccentricity. The non-zero 
eccentricity of the transiting planet may also betra y the pres- 
ence of an outer perturber in the system (see e.g. iFabrvckvl 
2009). Tidal interactions with the star are expected to cir- 
cularise the orbit of the transiting planet (|Zahn|[T977l ) but 
the presence of the second pla net in the wide r orbit may 
significantly delay this process (|Mardhndl2o"o71 ). 

The eccentricity of WASP-3b w as reported to b e sta - 
tistically indistinguishable from zero. iPollacco et al.l (120081 ) 
obtained e = 0.05 ± 0.05 and ISimpson et alj (|20ld )~~got 
e = 0.07 ± 0.0 8. Authors of both p apers adopted e = 0.0 in 
their analysis. iTripathi et al.l(|2010l ) adopted a priori circular 
orbit. We reanalysed the jo ined radial velo c ity da ta sets from 
IPollacco et"ai] l|200Sl) and ITripathi et"aH (|2010l ) (only out- 
of-transit measurements) to test solutions with WASP-3b 
orbiting in an eccentric orbit. We us ed the Systemic Con- 
sole software |Meschiari et alj|2009l ). The zero eccentricity 
Keplerian fit gives a reduced \ 2 = 3.3 with rms — 15.4 m 
s _1 . If the eccentricity e is set as a free parameter the best- 
fitting model results in e = 0.05 ±0.04 with a lower reduced 
X 2 = 2.3 and rms = 13.9 m s . Our results show that the 
orbit of WASP-3b may be circular within 2a and the ec- 
centric orbit cannot be excluded. The remaining dispersion 
of the radial velocity is typical for fast rotating late-type F 
stars |Saar. Butler fc Marcvlll998l ). 

4.3 Inner perturber 

The mass of the perturbing planet in an inner 2:1 resonance, 
M p , depends on the mass of the outer, more massive tran- 
siting planet M t and its apparent eccentricity e as follows: 

Mp = ~k Mt (1) 

l|Anglada-Escude et al.1 [2010). In this case M p should be 
treated as a upper limit because it was obtained assuming 
circular orbits. The hypothetical inner, less massive planet 
is expected to produce perturbations which should be visible 
in TTV of WASP-3b. 

To check if an inner perturber in the WASP-3 system 
may reproduce observed cyclic TTV, we generated 2232 syn- 
thetic O — C diagrams f or orb its close to the 1:2 resonance. 
Following iRabus et alJ ([2009) we neglected the light-time 
effects and assumed that the plan etary orbits are coplanar. 
We used the Mercury package (|Chamberslll999l ) and the 
Bulirsch-Stoer algorithm to integrate the equations of mo- 
tion for this three-body problem. Parameters of the transit- 
ing pl anet and its host sta r were taken from ITripathi et al.l 
(201Ch and IPollacco et al.1 (|2008h . respectively. The initial 
circular orbits were assumed. The semi-major axis of the 
perturber a p was between 0.0185 and 0.0215 au (±0.0015 
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Table 3. Outcr-pcrturbcr solutions which reproduce the ob- 
served O — C variation. a p denotes the semi-major axis of the 
perturbing planet, M p is its mass, P p is its orbital period, K p is 
the expected semi-amplitude of the radial-velocity variation and 
X? e d ' s *' le l° wes t value of reduced chi-square for direct model 
fitting. 



dp 


M p 








(au) 




(d) 


(m s ) 




0.0441 - 0.0443 


6-10 


3.03 - 3.05 


2.3 - 3.8 


2.4 


0.0493 - 0.0498 


10 


3.58 - 3.64 


3.7 


2.8 


0.0506 - 0.0511 


15 


3.72 - 3.78 


5.3 


1.5 



au away from the 2:1 resonance) in steps of 0.0001 au. The 
mass of the second planet M p was varied from 5 to 30 
in steps of 5 Mq^. The upper limit was calculated accord- 
ing to equation ([1} assuming e = 0.05 ± 0.04. The initial 
phase shift toward the transiting planet (in fact, an initial 
mean anomaly of the perturber), <j>, was between 0° and 
330° in steps of 30°. Computations covered 640 periods of 
WASP-3b, i.e. the time span covered by observations. Syn- 
thetic O — C diagrams were searched for periodicity with the 
ANOVA method and the mean amplitude of variation was 
determined. The procedure was repeated for the 3:2 reso- 
nance. 

Periodicity close to Pttv was not found near either reso- 
nance. That allowed us to eliminate the scenario where the 
inner perturber generates the observed TTV signal. 

4.4 Outer perturber 

Applying methods described in section 4.3, we generated 
19680 synthetic O—C diagrams assuming the existence of an 
outer perturber in the system. Both planets were assumed 
to orbit the star in initial circular orbits. We investigated 
cases with the second planet close to 2:1, 3:1, 3:2, and 5:3 
MMRs, ±0.0020 au away from each resonance. The mass of 
the perturber M p was varied from 5 to 50 Mq in steps of 
5 Mm. Solutions close to observed cyclic TTV were studied 
in detail by varying M p and <f) in steps AM P = 1 Mj and 
A(j> = 15°, respectively. 

We found a family of solutions close to 5:3 MMR re- 
covering the period close to P ttv for a p between 0.0441 and 
0.0443 au and M p between 7 and 10 Mfg. The exemplary 
synthetic O — C diagram is plotted in Fig. [4^,. We also found 
2 groups of solutions close to the 2:1 MMR. In the first one 
the perturber has M p of 10 Mj and a p between 0.0493 and 
0.0498 au. Fig.[4p shows the shape of expected signal which, 
despite being periodic, clearly shows a deformation of the 
amplitude. The second group of solutions may be found for 
a p between 0.0506 and 0.0511 au and M p of 15 Mq. As it 
is shown in Fig. |4j;, periodic variation may be modulated 
by a long-term periodicity significantly affecting the ampli- 
tude of variations. This long-term cyclic variability is caused 
by the resonant oscillations of the planetary orbits. Table [3] 
summarises parameters of proposed models. 

No solutions close to observed TTV variation were 
found near resonances 3:1 nor 3:2. Orbits far away from 
MMRs were analyse d using the ptmet code based on the 
pertu rbation theory (|Nesvornv fe MoTb idclli 2008; Nc svornvl 
l2009f ). A perturbing planet was put in initial circular or- 



bits with semi-major axis a p between 0.0370 and 0.1200 
au in steps of 0.0002 au. As the amplitude of TTV signal 
was found to scale nearly line arly with the perturber's mass 
IjNesvornv fe Morbidellil 120080 . M p was fixed and set equal 
50 Mgj . An initial phase shift toward the transiting planet <f> 
was between 0° and 315° in steps A<f> = 45°. A total of 3320 
systems were analysed but no configuration reproducing the 
observed variation was found. 

To quantify the quality of solutions collected in Table [31 
the synthetic O — C diagrams were directly fitted to data 
points. The chi-square test favours the third group of solu- 
tions where the smallest reduced Xrsd = 1-5 was obtained 
for a p = 0.0507. This value is significantly smaller than 
Xrod — 3.9 for a model without a perturbing planet. The 
O — C diagram with the best-fitting model is presented in 
Fig. [S] where the residuals are also plotted. The expected 
semi-amplitude of the radial velocity variation caused by 
the perturber is expected to be ~5 m s~* and thus unde- 
tectable in the published radial velocity measurements. The 
most probable model was found to be stable dynamically. A 
simulation was run for 10 6 yr and no significant changes in 
the eccentricities were detected. The eccentricities of tran- 
siting and perturbing planet were found to be below 0.005 
and 0.033, respectively. 

Such a WASP-3 system would not be completely un- 
usual with t wo planets close to the 2:1 MMR (see e.g. Gliese 
876 b and c lMarcv et"afll200lb but the first one with such 
close-in planets and the less massive planet being in an outer 
orbit. The origin of such a system would be interesting in a 
context of the large-scale orbital migration in young plan- 
etary systems. A Super-Earth (10 Mgj) is expected to be 
caught in a trap when it is very close to the edge of the wide 
gap opened in a disk by a Jupiter-mass gas giant. This trap- 
ping prevents the outer planet from being locked in the 2:1 
MMR. If a planet increases its mass by accretion, it may be 
released from such a trap and form a resonant configuration 
(IPodlewska fc Szuszkiewic3l2009h . This scenario cannot be 
rejected in the case of the outer perturber in the WASP-3 
system because the nowadays planet might be a rocky rem- 
nant core of a planet which has greater mass in the past - 
a hot Jupiter whose gaseous envelope was co mpletely evap- 
orate d due to the proximity to the host star l|jackson et al.l 
l2010h . 



5 CONCLUSIONS 

In this paper we show that transit timing of WASP-3b can- 
not be explained by a constant period of the exoplanet. We 
emphasise that this finding is based on only a few data 
points. However, one may put forward provisional hypoth- 
esis which assumes the existence of a perturber in the sys- 
tem. As a result of simplified 3-body simulations, two groups 
of configurations which reproduce the observed TTV were 
identified close to the outer 2:1 MMR and one group close to 
the 5:3 MMR. A model with the hypothetical second planet 
of the mass of ~15 Mgj and orbital semi-major axis of 0.0507 
au was found to be most likely. In this scenario both planets 
stay in the stable orbital resonance with a period ratio of 
2.02, i.e. very close to the 2:1 MMR. 

The expected radial- velocity semi- amplitude of the hy- 
pothetical second planet would be much smaller than the 
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Figure 4. Examples of synthetic O — C diagrams which reproduce 
the observed variations, generated for an outer perturber. Left 
panels - the unfolded TTV signal, right panels - the TTV signal 
folded to period close to Pttv 
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or disprove the claimed variation in the O — C diagram. 
A long-time baseline of observations is needed to charac- 
terise possible resonant oscillations in the TTV signal. More 
accurate radial-velocity measurements are also required to 
gain deeper insight into orbital parameters of the transit- 
ing planet. Additionally, a transit of the predicted second 
planet could be observable with a large ground-based or 
space-based telescope. 
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Figure 5. The O — C diagram with the best-fitting model which 
reproduces the observed variations. The mass of the perturber 
is 15 Mm and the semi-major axis of its initially circular orbit 
is 0.0507 au. Symbols are the same as in Fig. [2] Residuals are 
plotted in the bottom panel. 

stellar jitter of WASP-3. This finding reduces chances of a 
direct detection of the perturber's signal by further radial- 
velocity follow-ups. However, we showed that one cannot 
exclude the non-zero eccentricity of WASP-3b. Such an ec- 
centric orbital solution may be an indirect fingerprint of the 
second planet in a wider orbit gravitationally pumping the 
eccentricity of the transiting planet. 

Assuming the hypothetical second planet is also a tran- 
siting exoplanet, it would cause periodic flux drops of the 
host star. The depth of these transits would be in the 
range of 0.03-0.35 per cent (or 0.3-3.8 mmag) depend- 
ing on the adopted mean planetary density - for a rocky 
plane t (e.g. CoRoT-7b, iLeger et al.1 120091 ; lOueloz et atl 
2009) or an ultra-low -density hot Neptune (e.g. WASP-17b, 
Anderson et al 2010:. respectively. If one considers the mass 
of the hypothetical second planet and its proximity to the 
host star, the latter scenario seems to be less probable, caus- 
ing transits to be shallow and difficult to detect. 

As the TTV method requires many high-quality light 
curves, more data will be gathered for WASP-3b to confirm 
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